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Abstract 

We give a more complete calculation of 6 — > S7 decay, including leading log QCD cor- 
rections from mtop to M\y in addition to corrections from M\y to rrih- We have included the 
full set of dimension-6 operators and corrected numerical mistakes of anomalous dimensions 



in a previous paper [ 14 1. Comparing with the calculations without QCD running from mtop 



to Mvi/|12|, the inclusive decay rate is found to be enhanced. At mt = ISOGeV, it results 
in 12% enhancement, and for mt = 250GeV, 15% is found. The total QCD effect makes an 
enhanced factor of 4.2 at mt = 150GeV, and 3.2 for mt = 250GeV. 

PACS numbers: 12.38. Bx, 13.40.Hq, 13.20. Jf 
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1 Introduction 



Recently the CLEO collaboration has observed |1|] the exclusive decay B — K*'~f with a branching 
fraction of (4.5 ± 1.5 ± 0.9) x 10^^. A new upper limit on the inclusive 6 — 57 process is 
also obtained as B{b —>■ 57) < 5.4 x 10"'^ at 95% C.L.0. This has been a subject of many 
papers ^ ^, |^ recently. It has been argued that this experiment provides more information 
about restrictions on the Standard Model, Supersymmetry, Technicolor etc. Their results are 
found to be sensitive to the theoretical calculation of 6 — 57 decay. In order to reduce the 
theoretical uncertainty, a more accurate calculation of this decay rate is needed. 



The radiative b quark decay has already been calculated in several papers |P]-|TJ]. It is found 



to be strongly QCD-enhanced (e.g. a factor of 7 for nit = 80 GeV and Aqcd = 300 MeV in 



ref.[12 ). In other words, the strong interaction plays an important role in this decay. However, 



there are still some uncertainties in these papers. In ref . P| , the anomalous dimension matrix was 



truncated and the estimated uncertainty due to this truncation is less than 15%. Refs.[jTO|, |Tl 
did not contain a full leading logarithmic analysis either, although they included some of the 
terms neglected in ref. Ref. [jr^ is believed to be a more accurate result. But it did not 



include the QCD running from rritop to Mw- Since the top quark is found to be heavier than W 
boson( nitop = 174±10l;j2 GeV. [|I3]), a detailed calculation of this effect is needed . Ref.[Q does 
include this running, however there are some errors in the calculation of anomalous dimensions, 
which can lead to some changes in the final result. 

In the present paper, we recalculate the b —>■ s'j decay including QCD running from rritop 
to Mw, and correct the errors in ref. [Q, i.e. its anomalous dimension matrix. Furthermore, we 
use untruncated anomalous dimensions of QCD running from to mi,. 



2 Matching at ^ = rrit 

In Minimal Standard Model, we first integrate out the top quark, generating an effective five- 
quark theory. By using the renormalization group equation, we run the effective field theory 
^ See discussions in ref. . 
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down to the W-scale, at which the weak bosons are removed. Finally we continue running the 
effective field theory down to b-quark scale to calculate the rate of radiative b decay. To maintain 



gauge invariance, we work in a background field gauge JT5 



The charged sector of Standard Model Lagrangian is 



^cc = ~^l^'^'^92 [ u c t 



1,V 



wi 



+ 



1 



fJ' ' 92 



V2 Mw 
+ h.c. 



u c t 



MjjV 



u c t 



VM 



D 



R. 



Where V represents the 3x3 unitary Kobayashi-Maskawa matrix, Mjj and denote the 
diagonalized quark mass matrices, the subscript L and R denote left-handed and right-handed 
quarks, respectively. 

We first integrate out the top quark, introducing dimension-6 operators to include effects of 



the absent top quark. In ref.|Tj], an approximation was made to keep only leading order terms 
of 6 = M^/mj, neglecting the next to leading order charged current to W boson. To include the 
full set of dimension-6 operators, we have to pick up five more operators involving W bosons. 
Our operators now make a complete basis of dimension-6 operators. Higher dimension operators 
are suppressed by factor p'^/mf, where p"^ characterizing the external momentum of b quark etc. 



jo^ ~ ml. The basis operators are: 



O] 



LR 



^LR 



'^LR 



Qlr^ 



^rribSLD Br, 



Pt = ^^^'''§^F,-s,Y^'D^h,, 
Rl = tfi'gl{D''<l)+)<P^SLlM, 

Wl = ^^^'gl{D'^W';_)W^^Ll^.lal.bL. 
Wl = tfi'glW+^W!^SL$bL, 

Wt = ifi'glWlW^MD^^i. + iMbL. (2) 
Where SLD^'jubL stands for 'siD^'-^^hL + {D^SL)'yubL and the covariant derivative is defined as 

with g3 denoting the QCD couphng constant. The tensor appearing in P^'^ assumes the 
following Lorenz structure, the index A ranging from 1 to 4: 

Then we can write down our effective Hamiltonian 

Ueff = 2V2GFVt^y:,Y.Ci{^l)0,{^l). (4) 

i 

The matching diagrams are displayed in Fig.l and Fig. 2. The diagrams involving W bosons 
are introduced in addition to Goldstone boson ones, So the number of diagrams considered now 
is 2-times bigger than that in ref . . After tedious calculation we have 



- ^^l' - [ (1-5)3 + (1-5)4 ^""^^ 

/_!_ 15+ 1752_ 7^3 10^3 _ 4^4 X 



(5)3 (i-sy 

^-6-^5^ + 5'' , 5-35^ + 2(53 



(1-5)3 (1-5)' 



Cp3 = 0, 

Cr3 = 0, 

= C^r3 = Cy^i = 0, 

Cwl = = S/gj. (5) 

Notice that we have included the terms of all orders of 5 = / mf as far as the dimension- 
6 operators are concerned. However, to really achieve higher accuracy in 5, the dimension-8 
operators should be considered in the matching of Mw scaleQ. These coefficients are all from the 
finite part integrations of electroweak loops. Terms like log(/i^/m^) always accompanied by the 
infinity 1/e vanishes here, because of our matching scale fi = rrit. They will be regenerated by 
renormalization group running of electroweak in the next section. 



3 Running from nit to Mw 

The renormalization group equation satisfied by the coefficient functions Ci{fi) is 

/^|^^^(/^) = E(7n..Q(/^)- (6) 

Where the anomalous dimension matrix 'jij is calculated in practice by requiring renormalization 
group equations for Green functions with insertions of composite operators to be satisfied order 
by order in perturbation theory. Let Tq^ denote a renormalized n-point IPI Green function 
^We thank Professor R. Barbieri for making this point exphcit to us. 
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with one insertion of operator Oj. Then the anomalous dimension 'jij characterizing the mixing 
of operator Oj into Oj is determined from the renormahzation group equation for Tq^ , 



(7) 



Here (3 = fi{d/dfi)g, 7^ = {fi/m){d/dfi)m and wj^xt stands for the wave-function anomalous 
dimensions arising from radiative corrections to the Green function's n external lines. 

After evaluating the loop diagrams, we find the following leading order weak mixing of 
operators, with the Q, R part agrees with ref.[ffl]. 
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These mixing are all between operators induced by tree-diagram and loop-diagram. The vanish- 
ing log(/i^/mj) terms in the last section are regenerated here by renormahzation group equation. 
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The QCD anomalous dimensions for each of the operators in our basis are 
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(10) 



Comparing with ref.[Q, except the W part, there are still some differences in the anomalous 
dimension matrix, which may lie in omitting a factor of 1/2 in ref.[0] in calculating Feynman 
diagram like Fig.3. After these changes, the whole matrix can be easily diagonalized, and all 
eigenvalues are real, which is required to maintain hermiticity of the effective Hamiltonian at all 
renormalization scales. While in ref-Jl^ it can not. In their case, some eigenvalues are complex. 
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The solution to eqn.(^) appears in obvious matrix notation as 



C(/i2 



exp 



dg—r-— 

93(m) p[g) 



After inserting anomalous dimension (|8|-[T0D, we can have the coefficients of operators at = My/. 



And some of these operators change a lot from ref . [jl4[ due to our improvements. For details, see 
next section. 



4 Matching at /i = Mw 

In order to continue running the basis operator coefficients down to lower scales, one must 
integrate out the weak gauge bosons and would-be Goldstone bosons at yU = My/ scale. The 
diagrams are displayed in Fig. 4. In these new matching conditions, one finds the following 
relations between coefficient functions just above and below /i = My/'- 



Col,{Myr) 


= Coi^iM^ 


), 




= Coi^iM^ 


), 




= Coi^iM^ 


), 




= Cp.i(M+) 


+ 2/9, 




= Cp.,.(M+) 


-7/9, 




= Cpi,3(M+) 


+ 2/9, 


CpiA{M^) 


= Cpi,4(M+) 


+ 1, 


Cp2{M^) 


= Op|(M+) 


-Ch.2(M+)-3/2 


Cpz{M^) 


= Cp3(M+), 




Cp.{M^) 


= Op4(M+) 


+ 9. 



In addition to these, there are new four-quark operators: 

O3 = {sLa7''ha)[{uLpli,UL(3) + ■■■ + (bLp'Jf.bLfs)], 



with coefScients 



Oq = {sLal^hLp)[{uRp'^^URa) + ... + (hRp'^fJjRa)], 



Ci{Mw) = Q, 2 = 1,3,4,5,6, C2{Mw) = l. 



(13) 



To simplify the calculation and compare with the previous results, equations of motion(EOM) |]T6| 
is used to reduce all the remaining two-quark operators to the gluon and photon magnetic mo- 
ment operators and 0\ji. The effective Hamiltonian then appears just below the W-scale 
as 



AC 



w) 
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4Gf 

72 
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LR 
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(14) 



For completeness, we first give the explicit expressions of the coefficient of operator 
and at /i = M^, 
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They are expressed by coefficients a.t fi = rrit and QCD coupling a^. So it is convenient to utilize 
these formula. 
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In the previous paper |jT^, the higher order terms of M^^/mf are included in this stage by 
hands in order to match the previous workp[|-[p!^ when — ^ M]v. Therefore it is unnatural. 
While in this work, we keep higher order terms of M^/m'f from the very beginning at = rrit 
scale. If the QCD correction is ignored (by setting as{mt) = as{Mw) in eqn. ([T5D , ([T6|) ), our 
results would reduce to the previous resultsj^, exactly where the top quark and W bosons 
are integrated out together. This is required by the correctness of the effective Hamiltonian, and 
it is also a consistent check. 



If we rewrite our operators O^,^, as O7, Og like ref. |T2 



O7 = {e/lQ^'')m,SLa^'hRF^ 
Os = {g/16n^)mhSLa^''T%RGl,. 



(17) 



then 



Cs{Mw) 



The obvious differences from QCD correction to Cii^Mw) and CsiMw) can easily be seen from 



Fig. 5 and Fig.6. In comparison to ref.[|lj], the enhancement of coefficient of operator O7 is 
almost the same size, but the values for Og are quite different. Here the effect to Og is an 
enhancement rather than a suppression as in ref.|T^. At rrit = 150GeV, it is enhanced a factor 
of 40% in comparison to ref . |T3 . These changes come mainly from the corrections of anomalous 
dimensions described earlier. Since CriMw) and Cs{Mw) are both the input of the following 
QCD running from to m^. It will be is expected to change the final result. 



5 The B Xgj decay rate 

The running of the coefficients of operators from /i = Mp^/ to fi = ui}, was well described in 



ref. [12 . After this running we have the coefficients of operators at /i = scale. Table 1 
gives the numerical values of coefficients of operators O7, Og with different input of top quark 
mass. Here we use My/ = 80.22GeV, rrib = 4.9GeV. Both C-jinib) and C^inib) are enhanced in 
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comparison to values obtained by Misiak [|1^. The obtained values of Ce at the scale are 
-0.030, -0.036, -0.041 for A{fc-D = 100, 200 and 300 MeV, respectively. The corresponding values 
of C5 are 0.007, 0.008 and 0.009. They are just the same as ref. [|r^ . 

The leading order 6 ^ 57 matrix element of H^ff is given by the sum of operators O5, Oq 
and O7, 

< H,ff >= -2V2GFV*Vtb{CM + Qd[C,{fi)+3Ce{fi)]} < \0r\ > . (19) 
Therefore, the sought amplitude will be proportional to the squared modulus of 

C7^^^K) = C^{m,) + Qd [C^{m,) + 3^6 (m,)] (20) 

instead of \Cj{Mi,)\'^ itself. 



Following ref.[g-[|T2 



BR{B Xsj)/BR{B X^eV) ^ T{h s-f)/T{b ceV). (21) 



Then applying eqs.(p!9|),(|20|), one finds 



^ ^'-1' - ^^|C?"(rn.)P { 1 - ^-^fMn.^ " . (22) 

BK{B — > Xceu) TTg[mc/mh) \ J 

where g{mc/mf)) ~ 0.45 and f{mc/mi,) — 2.4 corresponding to the phase space factor and the 
one-loop QCD correction to the semileptonic decay, respectively ||1^. The electromagnetic fine 
structure constant evaluated at the b quark scale takes value as aqED^'^nb) = 1/132.7. Afterwards 
one obtains the B — > Xg'j decay rate normalized to the quite well established semileptonic decay 
rate. The results are summarized in Fig. 7 as functions of the top quark mass: The two upper 
solid lines represent the QCD-corrected ratio of the decay rates, corresponding to Aqce, = 100 
MeV and Aq^^ = 300 MeV, respectively. While the dashed lines correspond to the results 



obtained by Misiak|T2|. The QCD-uncorrected values are also shown. 



In this figure, we can easily see that, at rrit = ISOGeV, it results in 12% enhancement from 



Misiak's result|12], and for rrit = 250GeV, 15% is found. 
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6 Conclusion 



As a conclusion, we have given the full leading log QCD corrections (include QCD running from 
rritop to Mw), with whole anomalous dimension matrix untruncated. Comparison to the previous 
calculation three points are improved: 

(1) We have included the full set of dimension-6 operators. 

(2) We correct errors of anomalous dimensions in ref . . 

(3) We use untruncated anomalous dimensions in QCD running from to mh instead of 
truncated ones. 

In fact, point(2) makes an enhancement while point(3) leads to a suppression. For point(l) 
there is no significance change in final result. The total result does not change a lot, e.g. a 
suppression of 4% comparing ref. |T^ . Table 2 gives the results by different authors. 

The whole QCD-enhancement of the BR{B Xsj) makes a factor of 4.4 at rrit = 145GeV, 
and 3.2 at mt = 250GeV, when Aqcd = 200MeV. 

The branching ratio BR{B Xs'y)/BR{B — > XceV) ranges from about 3 x 10~^ to 6 x 10~^ 
as mt varies from 90GeV to 250GeV{AQCD = 200MeV). Although this result is not quite different 
from the previous calculations, our improvements lie in reducing some theoretical uncertainties. 



This improvements are important, since the anomalous dimension matrix of ref.0] is not a 
hermitian one, and its results are not faithful. 



The gluon magnetic moment operator Og is also enhanced from Misiak's result [|T2|, which 
is not the present paper's interest. 
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Table 1: Numerical results for coefficients of operators Ci{mb). 



mtop(GeV) 


cvirrib) 


csirrib) 




A=100 A=200 A=300 


A=100 A=200 A=300 


100 


-0.297 -0.324 -0.345 


-0.151 -0.163 -0.173 


110 


-0.307 -0.333 -0.353 


-0.155 -0.168 -0.177 


120 


-0.315 -0.311 -0.361 


-0.160 -0.172 -0.181 


130 


-0.323 -0.349 -0.368 


-0.163 -0.175 -0.184 


140 


-0.330 -0.356 -0.375 


-0.166 -0.178 -0.187 


150 


-0.337 -0.362 -0.381 


-0.169 -0.181 -0.190 


160 


-0.343 -0.368 -0.387 


-0.172 -0.184 -0.193 


170 


-0.349 -0.374 -0.393 


-0.174 -0.186 -0.195 


180 


-0.354 -0.379 -0.398 


-0.176 -0.188 -0.197 


190 


-0.360 -0.384 -0.402 


-0.178 -0.190 -0.199 


200 


-0.364 -0.388 -0.407 


-0.180 -0.191 -0.200 


210 


-0.369 -0.393 -0.411 


-0.181 -0.193 -0.202 


220 


-0.373 -0.397 -0.415 


-0.183 -0.194 -0.203 


230 


-0.377 -0.401 -0.419 


-0.184 -0.196 -0.204 


240 


-0.381 -0.404 -0.422 


-0.185 -0.197 -0.205 


250 


-0.384 -0.408 -0.426 


-0.187 -0.198 -0.207 



Table 2: Enhancement or suppression (-) of 6 ^ S7 decay rates obtained by Cho et al.[14], 
Misiak[12] and the present improved ones relative to earher result by Grinstein et al.[9] 





Grinstein et al. 


Cho et al. 


Misiak 


This paper 


QCD running Mw — > 


truncated 


truncated 


untruncated 


untruncated 


QCD running rrit Mw 


neglected 


performed 


neglected 


performed 


mtop=150GeV 


1 


+8% 


-6% 


+5% 


mtop=250GeV 


1 


+14% 


-5% 


+10% 
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Figure Captions 



Fig.l Leading order matching conditions at the top quark scale for the IPI Green functions in 
the full Standard Model and in the intermediate effective field theory. 

Fig. 2 One loop matching conditions at the top quark scale for the IPI Green functions in the 
full Standard Model and in the intermediate effective field theory. 

Fig. 3 One of the Feynman diagram in calculating Anomalous dimensions, with the heavy dot 
denoting high dimension operator. 

Fig.4 Matching conditions a,t fi — Mw for four quarks and two quarks IPI Green functions in 
the intermediate effective field theory and effective field theory below W scale. 

Fig. 5 The photon magnetic moment operator's coefficient Cr{Mw) for different top quark mass. 
The ones with and without QCD corrections are indicated by solid and dashed lines respectively. 
(A = 300MeV is used) 

Fig. 6 The gluon magnetic moment operator's coefficient Cs{Mw) for different top quark mass. 
The ones with and without QCD corrections are indicated by solid and dashed lines respectively. 
(A = SOOMeV is used) 

Fig. 7 BR(i? Xs'~f) normalized to BK{B XcCV), as function of top quark mass. The upper 
solid lines indicated our results for a full QCD correction. Dashed lines correspond to Misiak's 
results without QCD running from rritop to Mw- 



15 



This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9408351v2 



This figure "figl-2.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9408351v2 
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